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We present an efficient diagrammatic method to describe nonlocal correlation eflFects in lattice 
fermion Hubbard-like models, which is based on a change of variables in the Grassmann path 
integrals. The new fermions are dual to the original ones and correspond to weakly interacting 
quasiparticles in the case of strong local correlations in the Hubbard model. The method starts 
with dynamical mean-field theory as a zeroth-order approximation and includes non-local effects in 
a perturbative way. In contrast to cluster approaches, this method utilizes an exact transition to a 
dual set of variables. It therefore becomes possible to treat vertices of an effective single-impurity 
problem as small parameters. This provides a very efficient interpolation between band-like weak- 
coupling and atomic limits. The method is illustrated on the two-dimensional Hubbard model. The 
antiferromagnetic pseudogap, Fermi-arc formations, and non-Fermi-liquid effects due to the van 
Hove singularity are correctly reproduced by the lowest-order diagrams. Extremum properties of 
the dual fermion approach are discussed in terms of the Feynman variational principle. 
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I. INTRODUCTION 

One of the most successful theories of strongly corre- 
lated fermions on a lattice is dynamical mean-field the- 
ory (DMFT)^. Physically, this approach treats the lo- 
cal spin and orbital fluctuations of the correlated elec- 
trons in a correct self-consistent way , while the spatial 
intersite correlations on the lattice are neglected. The 
non-perturbative DMFT approach is successful, because 
a number of the most important correlation effects are in- 
deed related to local fluctuations. For example, DMFT 
describes correctly such phenomena, as the local moment 
formation in itinerant magnets^, some aspects of Kondo 
physics^, and the Mott insulator-to-metal transition on 
a lattice with a large connectivity in high-dimensional 
materials^. 

On the other hand, there is increasing evidence that 
the non-locality of spatial correlations plays an impor- 
tant role, particularly for the Luttinger liquid physics 
of low-dimensional correlated systems^, d-wave pairing 
in quasi two-dimensional cuprates^, and non Fermi- 
liquid behavior due to van-Hove singularities in two- 
dimensional systemaii^i^. Moreover, angle-resolved pho- 
toemission spectra of three-dimensional ferromagnetic 
iron shows appreciable k-dependent self-energy effecta^^. 

The most obvious generalizations of DMFT that takes 
into account the short-range non-local fluctuations are 
the so-called cluster DMFT approximations, in real or k- 
spaceiiii^. In these methods, correlations are assumed 
to be localized within a cluster including several lat- 
tice sites. Cluster methods do catch the basic physics 
of c?-wave pairing and anti-ferromagnetism in high-Tc 
superconductorsi^ii^ and the effects of inter-site Coulomb 
interaction in transition-metal oxidesi^. At the same 
time, the complicated k-dependence of the self-energy 



close to the Fermi surface, giving rise to Luttinger liquid 
formation is related to long-range fluctuations and there- 
fore cannot be described within cluster approaches. For 
the same reason, cluster methods hardly can handle the 
effects due to van-Hove singularities or nestingii^. An- 
other drawback of the cluster methods is that the specific 
choice of the cluster and corresponding self-consistency 
condition is not unique. Different self-consistency con- 
ditions (e.g. DCAi^ and free-cluster CDMFT") or 
periodization schemes (e.g. self-energy and cumulant 
periodizatioi*i^) can result in physically different solu- 
tions. For example, the critical temperature of the d- 
wave superconducting transition of the doped Hubbard 
model is different in DCA calculations^^ and for then 
2x2 free cluster 

The present paper is devoted to an alternative exten- 
sion of DMFT, which operates with a single-site impurity 
problem and treats spatial nonlocality in a diagrammatic 
way. 

Let us first recall the key DMFT equations. Formally, 
the assumption of local correlations means that the en- 
vironment of a correlated atom can be replaced with a 
Gaussian effective medium. Consequently, the lattice 
problem reduces to the impurity problem. The later is 
described by the effective impurity action 



A, 



(1) 



where Sat is an action of the isolated or bare atom, and 
the second term is the hybridization due to the rest of 
the lattice. An important property of the DMFT ap- 
proach is that this hybridization function has non-trivial 
frequency dependence, so that the approximation catches 
the physics of local fluctuations of spin, charge, and or- 
bital degrees of freedom. For example, it is vital, for the 
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description of Kondo physics-. 

It is obvious that the impurity problem is much simpler 
than the original lattice one. Nowadays, a number of nu- 
merically efficient impurity solvers are available. In par- 
ticular, these solvers allow one to calculate the Green's 
function of the impurity problem g^^^a on the Matsub- 
ara frequencies axis. This is the only property of the 
impurity problem entering in the DMFT self-consistent 
equations. The DMFT approximation for the Green's 
function of the initial lattice problem corresponds to the 
following expression 

r<DMFT _ 1 fcy\ 

One can see from this equation that the self-energy is 
local in DMFT, since the momentum dependence of is 
not renormalized. The hybridization function A satisfies 
the self-consistentency condition of DMFT, 



where Gr=Q — N^^ Gk is the local part of the Green's 
function ^ of the lattice with N sites. 

In order to understand the main idea of the present 
work, let us first describe in a simple way the DMFT 
condition fSj. If we consider the case of a truly Gaussian 
system then the DMFT approach becomes exact. For 
this case, equation ((S]) is trivial. Indeed, to obtain the 
impurity problem for the site j, one integrates out truly 
Gaussian degrees of freedom for other sites. This exact 
procedure does not change the properties of the electron 
motion at the site j, so the local part of the Green's func- 
tion before integration must equal the Green's function 
after the integration, Gr=q = g. Turning back to the 
general case of a non-Gaussian ensemble, we note that 
among different properties of the impurity model, the 
DMFT scheme uses only the local Greens's function g^^^- 
Once g^^ is known, the approximation does not differ 
between Gaussian and non-Gaussian cases. Therefore, if 
a certain equation for g^^ is established for the Gaussian 
limit, it must also remain valid for the general case. 

As it follows from the previous discussion, the DMFT 
equations are essentially the formulae for the Gaussian 
limit, renormalized in terms of the Green's function of the 
impurity problem. It turns out that the resulting theory 
works well, not only in the case of weakly interacting 
systems, but also in the atomic limit case, which is very 
different from a Gaussian system. A good interpolation 
between the two different limits is a key advantage of the 
DMFT approach. 

Starting with the above interpretation of DMFT, it 
is natural to discuss a possible extension of this theory. 
Such an extension should be based on the perturbation 
series near the Gaussian limit, renormalized in terms of 
the impurity problem. The lowest-order term of such a 
theory should restore the DMFT result, whereas higher- 
order corrections would describe spatial non-locality. A 
properly constructed theory of this kind would describe 



both short- and long-range fluctuations and will not suf- 
fer from the periodization problems of cluster DMFT. 

Unfortunately, the straightforward construction of 
such an extension meets serious difficulties. The prob- 
lem is that the extension is not unique. Beyond DMFT, 
there are many ways to choose the renormalization pro- 
cedure, to define the hybridization function for the impu- 
rity problem and other quantities. One can formulate the 
major requirements for the desirable non-local correlated 
theory, they include: 

• at least in the Gaussian and atomic limits, the the- 
ory should become a regular series around DMFT, 
with an explicit small parameter; 

• the basic conservation laws should be fulfilled in 
the theory; 

• the choice of hybridization function should be op- 
timal, in a certain sense; 

• there should be good practical convergence of the 
series: the leading corrections should capture most 
of the non-local physics; 

• last but not least, the equations of the theory must 
be easy enough for practical calculations. 

There have been several previous attempts to construct 
a proper theory of this kind^^'^^'''^*. These approaches re- 
quire a solution of ladder-like integral equations for the 
complete vertex F and the subsequent use of the Bethe- 
Salpeter equation to obtain the Green's functions. The 
first step exploits the vertex part of the effective impu- 
rity problem, whereas the second step uses just the bare 
interaction parameter U. We do not know of detailed 
tests of these approache o^^'^^i^^ , but we suspect that the 
presence of bare U in the theory makes it suitable for 
the metallic phases only. We also note that ladder-like 
integral equations are hard for practical calculations. 

In this paper, we describe in detail a formalism fulfill- 
ing all the criteria from the above list. A preliminary 
version of this method was published in Refii^. The 
method is based on the transition to the new set of vari- 
ables, called the dual ensemble. The procedure utilizes 
a Hubbard-Stratonovich transformation for the Gaussian 
part of the action. Several years ago, this trick was first 
proposed for classical fluctuation fields^. For a strong 
coupling expansion of the Hubbard model around the 
atomic limit without hybridisation function, the equiva- 
lent Hubbard-Stratonovich transformation has been pro- 
posed in different papers^^'^^. A similar procedure for 
fermions with general non-local interactions have been 
discussed recently^^. Also we would like to mention a 
much earlier work^l for classical fields. Although, it used 
a different formalisro^lj the resulting diagram series re- 
sembles ours. 

The paper is organized as follows. Section II is de- 
voted to the general theoretical framework. Section III 
describes the application of the non-local theory to the 
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problem of the antiferromagnetic pseudogap and the for- 
mation of Fermi arcs in the two-dimensional Hubbard 
model for high-temperature superconducting cuprates. 
In the Appendix A we discuss how the many-particle 
excitations for the initial and dual system are related. In 
the Appendix B the functional minimization derivation 
of the self-consistent DMFT condition is discussed. 



II. DUAL FERMION FORMALISM: BEYOND 
DMFT 

A. Definitions 

We start from the two dimensional Hubbard model 
with the corresponding imaginary-time action 

C*] = ^ (e/c - - i^) clkaCu;ka + U ^ / m-^rnilrdT. 

(4) 

Here f3 and fj, are the inverse temperature and chemical 
potential, respectively, oj = {2j + l)7r//3,j = 0, ±1,... 
are the Matsubara frequencies, r is imaginary time, 
a =t, I is the spin projection. The bare dispersion law is 
ek — — 2t(cos kx + cos ky), c*, c are Grassmanni variables, 
n-ia-T ~ c*^T'^ia-Ti where the indices i and k label sites and 
quasi- momenta. 

In the spirit of the DMFT, we introduce a single-site 
reference system (an effective impurity model) with the 



action 

/■/? 

Simp ^y^i^ui- fJ.-i(^)clj „C^,a + U / n^TnirdT, (5) 

where is a yet undefined hybridization function de- 
scribing the interaction of the effective impurity with 
a bath. We assume that all properties of the impurity 
problem such as single-particle Green's function g^, and 
higher momenta can be calculated. In particular, we will 
use the forth-order vertex 71234 — 5rr522^(Xi'2'3'4' — 
gvv92'z' + 5i'3'ff2'4')53^35^4 (l^^re, x is a two-particle 
Green's function of the impurity problem, and indices 
stand for a combination of a and uj, for example gut 
means gat.ujt.a-,, ,u>-,,)- 

Our goal is to express the Green's function G^k and 
other properties of the lattice problem of Eq. ^ via the 
quantities for the impurity problem. 



B. Dual variables: exact formulas 

Since A is independent of fc, the lattice action (|4]) can 
be represented in the following form 

S[C, C*] = ^ Sirnp[Ci,C*] - ^(A^ - efc )c* fe^C^^fe,^ . (6) 
i ujkcT 
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We utilize a dual transformation to a set of new Grassniann variables /, /* . The following identity 

is valid for arbitrary complex numbers A and a. We chose = (A^^ — e^) for each set of indices fc, a. The quantity 
a remains yet unspecified, but we require it to be dispersionless; a — a^^^a- 

With this identity, the partition function of the lattice problem Z — j e^'^'^^^ l2?c*2?c can be presented in the form 
^ = 11 e-^^"'"' '^'^'^Vf*VfVc*Vc, where 

S[c, c*, /, /*] = - E.fe In («L(A. - efc)) -I- ^™p[c., c*] + 

(8) 

As a next step, we establish an exact relation between the Green's function of the initial system Gr-r' = — < 
TcriC*,^, > and that of the dual system Gf^^, = — < Tfrif*,^, >■ To this aim, we can replace Ck — > £k + Se^jk 
with a differentiation of the partition function with respect to Se^jk- Since we have two expressions for the action ^ 
and ([8]), one obtains 

G^.k = (A^ - ek)-^a^,Gtj'a^,{A^ - e,.)"! + (A^ - ek)-\ (9) 

Similar relations hold also for higher-order momenta, as Appendix A describes. 

The crucial point is that the integration over the initial variables c* , Ci can be performed separately for each lattice 
site, since a is local and J^k (/fe^fc + c^/fc) = J2i + c*/i). For a given site i, one should integrate out c*, Cj from 
the action that equals 



(10) 
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We finally obtain an action S depending on the new variables /, /* only; 

where zl'^^ — J e~'^^"^pt^^ '^^^Vc^Vci, and the dual potential Vi = t^[/*, /i] is defined from the expression 



The Taylor series for V[fi, /*] can be obtained from the 
expansion of this definition in powers of /i, /*. One can 
see that (|12p defines V in such a way that this series 
starts from the quartic term, oc f*f*ff. Later on we 
take, for convenience: 

au;<T = gj\ (13) 
as it gives a particularly simple form of V. In this case the 
leading term in V is — 171234/1 /^/sA- Further Taylor 
series terms yield similar combinations including 7("'' of 
higher orders. 

Thus we see that in the dual action, the interaction 
terms remain localized in space, but are they non-local 
in imaginary time, since, for example, 7^'*) depends on 
the three independent Matsubara frequencies. Except 
for this point, the action (fTT|) formally resembles (j4|). 

There is a point which is worthwhile to discuss here: 
one can formally apply the transformation ([7]) with some 
new hybridization function to the dual system and 
thus obtain a sequence of changes to new variables. It is 
useless, however, since mathematically, these transforma- 
tions form a group. It is easy to show that any sequence 
of the transformations ([7|) corresponds to a single change 
of variables with a certain A. Moreover, there is an in- 
verse change of variables, that allows to obtain S[c, c*] 
back from the S[f,f*]. It is given just by Eq.d?]) with A 
replaced with aA^^. 

C. Gaussian approximation for dual ensemble 
and the relation to DMFT 

Since the transformation from the initial system (jj]) to 
the action (jlip contains no approximations, it is equally 
hard to describe exactly the properties of c* , c fermions 
as thereof /*, / dual-fermions. The main idea of switch- 
ing to the new variables is that, for a properly chosen 
A, correlation properties of the /*, / system are simpler 
than for the c*,c original model. In other words, the 
magnitude of the nonlinear part in the dual action can 
be effectively decreased by the proper choice of A. To il- 
lustrate this statement, let us just neglect V in (fTT|) . We 
denote the Green's function for such Gaussian approx- 
imation for the dual potential with calligraphic letters. 
The expression (fTTj) corresponds to 

Gif = ((Ac. - eu)-' + g^y^ 9^. (14) 



Being combined with the identity ([9]), this gives the for- 
mula 

^^,fc = (.gj'+A^-efc)"'. (15) 

One can recognize that this is exactly a DMFT expression 
for the Green's function. Therefore we conclude that for a 
properly chosen A already a Gaussian approximation for 
the dual potential yields a reasonable result, as DMFT 
does. It is important to point out that DMFT works 
well for the whole range of the parameters. In contrast, 
the Gaussian approximation for the atomic limit of the 
initial model ^ makes no sense. In that aspect, the dual 
potential V is indeed smaller than U . 

An argumentation can be presented to justify that the 
DMFT value of A is a proper choice for the Gaussian 
approximation [T5 |) . One of the reasons is described 
in the Appendix B. It turns out that Feynman minimiza- 
tion criterion for the Gaussian trial action, been formu- 
lated for the dual ensemble, gives exactly the DMFT hy- 
bridization function. Another argument is presented in 
the following subsection. 

Once the dual potential is taken into account, it yields 
a correction to the DMFT result. It is useful to introduce 
the dual self-energy 

^ dual ^ Q dual ^ '^dlal'^ (16) 

and the correction to the DMFT self-energy 

T.' = g-^-G-^. (17) 

With these quantities, we can reexpress the exact relation 
([9| in a particularly simple form 

We note that this expression relates quite different quan- 
tities: Srfua; and S' characterize the corresponding lat- 
tice problems and carry, in general, both momentum- and 
frequency-dependence, whereas g comes from the impu- 
rity model and is local in space. 

D. Diagram series: general properties and the 
choice of hybridization function 

The main idea of our method is to consider a diagram- 
matic expansion with respect to the dual potential V . 
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FIG. 1: Various diagrams for S''"'"'. Dia grams a, a', and a" 
are vanished by the condition 1191 

O 

FIG. 2: Two simple diagrams for Baym functional 
^duat[Gduat]- Fuuctioual differentiation of these diagrams 
with respect to Gduai produces diagrams a and b for self- 
energy. 

We will later demonstrate that already low-order dia- 
grams of such a series bring an important information 
about non-local correlations. The basic reasoning for this 
is presented in the previous subsection: since the value 
of V is in certain sense small, the first few terms of the 
perturbation series with respect to V can make sense. 
More detailed discussion about the small parameters of 
the theory are presented in the next sections; let us first 
present the general properties of the diagrams under con- 
sideration. 

The rules of diagram construction are quite similar to 
the usual Matsubara diagram technique. The only dif- 
ference from the standard perturbation scheme is that 
the interaction operator V is not purely of the 4-th order 
form f* f* f f , and therefore vertices in the diagrams are 
not necessarily four-leg, but may formally have any even 
number of legs. For the choice these vertices are 

essentially 7^"^. They are connected with the lines being 
the dual Green's functions. Some of the diagrams con- 
tributing yo the dual self-energy are presented in Figure 

We use the skeleton diagrams with renormalized 
Green's functions, so that the lines are complete GduaU 
and not Qduai- The reason to use the skeleton-diagram 
expansion for the dual self-energy is that it makes pos- 
sible to obtain conserving theories, similarly to conven- 
tional diagram technique^^. The Baym criterion of a con- 
servative theory is the existence of a functional of the 
Green function <1>[G] such that |^ = S. Once this func- 
tional is described by certain skeleton diagrams, taking 
the derivative means just cutting the lines in that dia- 
gram. For example, the diagrams (a) and (b) for the 
self-energy come from diagrams (a$) and (b$), shown in 
Figure([2]) (of course, care should be taken of the numeri- 
cal factors). Second-order differentiation with respect to 
G gives the two-particle quantities. Such a procedure au- 
tomatically produces a theory fulfilling the conservation 



laws for energy, momentum, particle numbers etc. 

In our consideration, the usage of skeleton diagrams 
describes a corresponding Baym functional ^duai[Gduai\ 
with the functional derivative being T^duai- Therefore, 
it produces a conservative approximation for the dual 
ensemble. Then it turns out that the exact transfor- 
mations © and (|A5p give a conserving description of 
the initial system. Simply, the conservation laws imply 
certain selection rules for G and F, and Q lASP clearly 
preserve those selection rules during the transformation 
from dual to initial quantities. More precisely, the con- 
serving character of an approximation in fact means that 
there exists some conserving dual action /*], exactly 
corresponding to this approximation. Since there is a 
one to one correspondence between S[J,f*] and ^[c, c*] 
(see the end of Section III Bp , we conclude that the initial 
system described by a certain S[c, c*] is also conserving. 

Until now, the hybridization function A was formally 
not specified. Now, we establish a condition for A that 
corresponds to a particular condition for the diagram- 
matic series. Let us again consider the DMFT. Suppose 
that we want to obtain the DMFT result without DMFT 
loops, that is using Ac^ not fulfilling (IBSp . Formally, it is 
possible: one should just sum up all the diagrams con- 
taining a single vertex (diagrams a, a', a" etc.). Since 
these diagrams give exactly the DMFT self-energy, such 
a procedure would indeed recover the DMFT result for 
an arbitrary hybridization function. The special DMFT 
choice of A just allows to eliminate such an infinite sum- 
mation, since (|B6p eliminates all the diagrams containing 
a simple closed loop. It is reasonable to keep this prop- 
erty in higher approximations, that is to require 

G'ZU = (19) 

as a condition for A. Then, all the diagrams with simple 
closed loops drop out from the calculation. Note that 
these diagrams however should be taken into account 
while taking the functional derivatives. For example, the 
DMFT vertex part Tduai = 7*''''' comes out from the dif- 
ferentiation of diagram (a). Finally, the condition (fT51) 
obviously passes into (jB6p at the DMFT limit. There- 
fore, until the corrections to DMFT are significant, one 
can approximate A„ with the DMFT hybridization func- 
tion. 

The vanishing of the closed loops seriously reduces the 
number of the low-order dual diagrams. In most of the 
practical calculations presented below we consider a sin- 
gle diagram (b). It is clear that any reasonable expan- 
sion starts from this perturbation, and that this diagram 
already incorporates some non-local physics. The corre- 
sponding formula for the dual self-energy reads (spin and 
orbital indices are omitted): 

Y^dual 1 \ ^ /.^ /^/^"^^ ^dual ^dual /-idual 

(20) 
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E. Causal properties 

Beyond conservation laws, the Green's function should 
be causal. The retarded Green's function G^{t), that is 
an analytical continuation of Gr to the real-time axis, 
should vanish for negative time: 

G^it < 0) = 0. (21) 

In the Fourier representation, condition ([2T|) implies 
the analyticity of G^ in the upper complex plane, as this 
follows directly from the definition of the Fourier trans- 
form. The inverse is also true. If the Fourier transform 
of a function is analytical in the upper-plane, the func- 
tion is causal. To prove this statement, it is enough to 
transform the integration contour of the inverse Fourier 
transform away from the real axis. 

Frequently, the causality principle is associated with 
the positiveness of the imaginary part of the Green's 
function in the real-frequency domain. For dual Green's 
function, this can lead to certain misunderstanding. It 
is clear from condition (jl9p that the imaginary part of 
Qduai (ja^jjj^Q^ ]-,g always-positive. However, this issue is 
purely formal. Condition (|2ip itself does not imply that 
ImG[^ is positive. A trivial counter-example is the func- 
tion — G^. It fulfills pTjl . and has an always- negative 
imaginary part. We will argue the same for G'^""'. It 
fulfills (HJ). Although, its imaginary part is not always- 
positive. 

Let us illustrate this statement at the zeroth order of 
the theory, single-site DMFT. It has been provertiiiiii^ 
that this theory is causal, so Q and g fulfill (pij) . One 
can easily check, from the expressions (|141 [TS]) . for the 
case of DMFT, a simple relationship holds; tj''™' = G — 
g. It is immediately clear from this formula, since G 
and g are causal, (^''""' also fulfills (pi]) . Note again, 
both condition ^ is fulfilled in DMFT, and -Img^""' 
is therefore essentially non-positive. 

Let us now consider the dual-fermion theory beyond 
DMFT. We will show that, if the hybridization function 
A is casual, the resulting Green's function is also causal. 
First of all, the casuality of A is inherited from g and 7^"^ 
Therefore, the dual system is characterized by the casual 
bare propagator g^-^'^^ = g — g and casual interaction 
operator. Therefore, the theory with skeleton diagrams 
results in a causal Qduai26^ Finally, it should just be 
proven that the casuality G^""' means the casuality of 
G. The later statement follows from the exact relation 
Indeed, since a = g^i does not have zeros in the 
upper-plane, g~^ is analytical. The same is true for the 
quantity (A — e)^^. Therefore, the entire right-hand side 
of ([9]) is analytical in the upper-plane. This implies the 
causality of G. 

In the calculation procedure described below, we al- 
ways start from a causal A and change it iteratively to 
deliver condition We will argue that such an iter- 

ation procedure preserve the causality of A. Therefore, 
the entire theory is causal. 



Finally, let us recall the issue of the positiveness of 
— ImG in the complex upper-plane. Actually, this is re- 
lated with the positivity of the residuals, as it follows 

from the Lehmann representation G = — . 

Here, the causality follows from the positivity of 0, 
whereas the requirement Z > ensures that — ImG > 0. 
For our theory, we were not able to prove the positivity of 
the residuals formally. However, we do not consider this 
as a serious drawback, since our practical calculations 
always produce undoubtly positive residuals. 

F. Small parameter in the extreme cases 

An important property of the DMFT approach is that 
it becomes exact for the two opposite cases of a non- 
interacting Gaussian system and of an extreme strong- 
coupling limit corresponding to the atomic limiti. The 
dual-fermion formalism inherits this property; moreover 
the corresponding smallness appears in the diagrams in 
a simple form. Let us first consider the strong-coupling 
limit Ek — *■ 0. It is useful to estimate the DMFT dual 
Green's function g, defined by the formula and the 
condition (fT9|) . For a pure atomic limit = 0, the 
Green's function is local, Gr^o = 0. However, the lo- 
cal part of the Green's function also vanishes due to the 
condition (HH). Formally, ^ ^ as A ^ 0. The small- 
ness of e and A allows the approximate estimation of 
the dual Green's function near the atomic limit. It gives 
Guik ~ guj^kgui- Since the DMFT is almost exact near 
the atomic limit, the same estimation is valid for G''""'. 
Consequently, near the atomic limit the lines in the dual 
diagrams carry a small factor ek- 

On the other hand, for the opposite weak-coupling 
limit [/ — > 0, the vertex parts of the impurity problem 
can be estimated as oc t/,7('5) cx U^, etc. There- 
fore, for the weak- coupling limit the vertices in the dual 
diagrams are manifestly small. 

The presence of a small parameter in these two limits 
does not guaranty a good interpolation between them. It 
should however be mentioned, that the scheme performs 
well if the corrections to DMFT are small: for this case 
we deal in fact with a perturbation series around DMFT. 
The validity of the method for more general situations 
should be checked in practical calculation. This practical 
validity depends on the particular choice of diagrammatic 
approximation for Y.(iuai- In this context, it is worth to 
discuss the choice of hybridization function A. 

G. Calculation procedure 

In practical calculations the solution was obtained it- 
eratively, similarly to the DMFT loop. The iterative 
scheme is presented in Figure |3l It includes the big 
(outer) and small (inner) loops. The small loop is de- 
voted to obtain the dual Green's function and self-energy, 
given the solution of the impurity model with certain 
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CT-QMC 



\^dual 



Gdual 



FIG. 3: (color online) The scheme of calculation. The calcu- 
lation includes "big" and "small" loop, marked with red and 
black lines, respectively. The small loop is to determine the 
renormalized dual Green's function G"*""' in a self-consistent 
way, for given A,g, and 7'"-'. The big loop is to determine A. 
Only the big loop requires a solution of the impurity problem. 



A. It starts from some guess for J^duai, for instance 

^flai = 0- The dual Green's function {Gduai~^duai)~^ is 
substituted in formula (|20p to produce a new estimation 
for T.duai- The procedure is repeated until converging 
results are reached. 

The big loop is very similar to the DMFT iterative 
procedure. We start with some initial guess for A and 
solve an impurity model. We use the weak-coupling CT- 
QMC solver—, which produces both the Green's function 
g and the 4-point vertex 7(4) in the frequency domain. 
Then we perform the inner loop to obtain Gduai (this 
step is not necessary in DMFT, since it uses the bare dual 
Green's function Qduai)- Finally, we take a new guess for 
the hybridization function 



A^ ^ A^ + C 



1 



(G*",lo)-i+.g: 



(22) 



and repeat the self-consistent procedure. A value of the 
parameter ^ < 1 was chosen to ensure better convergence. 
The last formula is organized in such a way that (i) its 
fixed point clearly satisfies the condition (|19p and (ii) for 
'^duai = it passes into the DMFT update formula A^^ 
Aij + £,{QZ\-=o ~ 9<^)- Of course, only the requirement (i) 
is actually necessary, so that formula is not unique. 
In particular, it is useful to consider an update 



(23) 



One can easily see that an update conserves causal 
properties of A, do that the convergence of the iteration 
process (^51 proves the causality of the result. Such a 
convergence indeed takes the place for the calculations 
presented below. Note also that near the fixed point 
Gf,"°lo = formula ^ passes into so that there is 
no much practical difference between these two formulas. 



III. APPLICATION TO THE HUBBARD 
MODEL 

In the next sections, we present the results of our cal- 
culations for the 2D Hubbard model. We start with the 
half-filled case with next-nearest neighbor hopping t' = 
lattice. We compare our data with direct QMC simula- 
tions on a finite Hubbard lattice, which are relatively 
simple due to the absence of the sign problem for the 
half-filled Hubbard model. 

Properties of the half-filled Hubbard model are well- 
known and are mostly related to the antiferromagnetic 
phenomenon and Mott metal-insulator transition. Lo- 
cal magnetic moment on atoms are formed and tend to 
ordered into an antiferromagnetic lattice due to the effec- 
tive super-exchange coupling. At zero temperature, the 
antiferromagnetism arises already at U = O"*", because 
of the perfect nesting. At finite temperature, the true 
antiferromagnetism is destroyed by the long-range fluc- 
tuations. However, short-range antiferromagnetic corre- 
lations are still present. Short-range antiferromagnetic 
ordering manifests itself as the strong pseudogap in the 
local electron spectral function. 

We consider the system with t = 0.25 at inverse tem- 
perature P = 20, with different values of U. Since the 
temperature is relatively high, it is enough to use the ref- 
erence data obtained just for the 8x8 lattice QMC sim- 
ulation, with subsequent maximum-entropy continuation 
of the data to obtain local density of states (DOS). The 
result for paramagnetic calculation is presented in Figs, 
m [5] (thin solid line). These results show that the narrow 
antiferromagnetic pseudogap is formed at approximately 
U = 1.0. For larger U, the DOS contains also a wider 
Mott gap, having a halfwidth of about U/2. At U ^ 2.0, 
the system shows essentially Mott-insulator DOS; the ef- 
fect of antiferromagnetism in this case consists in the 
sharp shoulders of the Mott gap. 

To understand better the physics of the half-filled Hub- 
bard model, it is worth to analyse the behavior of the 
electronic self-energy S. At small U, this is a small regu- 
lar correction to the dispersion law . It follows from the 
weak-coupling analysis that ImS is strongly anisotropic 
in this regime, with peaks near [0, ±7r], [±7r, 0] points. In 
contrast, for the truly antiferromagnetic gap, Sfc would 
have a pole at the Fermi surface. The residue of this pole 
is the same at all points of the Fermi surface. For large 
enough U this pole is somehow shifted from the real- 
frequency axis due to long-range thermal fluctuations, 
but the qualitative picture remains the same: a sharp 
peak in ImE, with almost constant magnitude along the 
Fermi surface. 

It is well-known that doping changes the physics of the 
Hubbard model substantially. First of all, already a few- 
percent doping suppresses the antiferromagnetism. At 
higher doping values there is a trend to d-wave supercon- 
ductivity. A superconducting phase has been obtained in 
various cluster-DMFT calculation o^^i^^ near the optimal 
doping of about 15%. This agrees well with the phase 
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diagram of high- Tc cuprates^. The pseudogap forma- 
tion in the doped Hubbard model was first analyzed by 
the cluster DMFT method (more specifically, Dynamical 
Cluster Approximation) in Ref. '291. For further appli- 
cations of the DCA to the 2D Hubbard model, see Refs. 
[l^l30H3lll32l . In the following consideration, we will not 
discuss the superconductivity itself, but we will address 
the so-called Fermi arc phenomenon. Essentially, this 
is an anisotropic destruction of the Fermi surface in the 
pseudogap regime. Only the parts of Fermi surface near 
the nodal direction remain well-defined at low tempera- 
ture. In the anti-nodal direction, the spectral function at 
the Fermi level is vanishingly small. 

A methodological difference between the doped and 
the undoped cases is that the sign problem makes di- 
rect lattice simulations away from half filling practically 
impossible^. Therefore, the reference point can only be 
the results of different approximate schemes or the ex- 
perimental data. 



A. Undoped case: translationally-invariant solution 

First, we discuss the result of the dual fermion investi- 
gation without a spontaneous symmetry breaking which 
means that the impurity problem is assumed to have no 
spin-polarization. The data presented in this chapter 
have been partly discussed previously as a Brief Report^^. 

The translationally-invariant DMFT predicts a Mott 
transition at rather high value U > 3.0 (for a bandwidth 
W = 8t = 2.0). It is important to point out that the 
density of states at the Fermi energy is independent of 
U within the entire Fermi-liquid phase. This is a conse- 
quence of the locality of the self-energy in DMFT. There- 
fore, for U ~ 1.5 — 3.0, the approximation predicts a 
three-peak DOS which consists on two Hubbard bands 
at ±C//2 and a Kondo-like central peak providing the 
'pinned' value of DOS at Fermi level. 

This behaviour is inconsistent with the reference data 
described above. Actually, those data do not show a 
three-peak structure, because of the antiferromagnetic 
pseudogap. Besides antiferromagnetism, the DMFT suf- 
ficiently overestimates the critical value of U for the Mott 
transition: according to the reference data, the the sys- 
tem shows DOS of the Mott-insulator nature already at 
U « 2.0 (see Figurei]). 

Let us take the leading dual diagram (b) into ac- 
count. The corresponding data is presented in Figs. IH 
[SJ Since the self-energy is not local anymore, there is no 
pinning at Fermi level, and the Kondo-like peak disap- 
pears. Furthermore, the self-energy momentum depen- 
dence agrees well with the qualitative picture described 
above. The upper panel of Figure [5] presents contour 
plots for ImE^=o,fc a.t U — 1.0 and U = 2.0 (the data are 
obtained by a polynomial extrapolation from the Mat- 
subara frequencies). The value of ImS^^^o./t grows dra- 
matically as U changes from 1.0 to 2.0. At larger U, there 
is an expected sharp non-Fermi liquid peak in ImE^^^o./t 
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FIG. 4: (color online) Local Green's function at Matsubara 
frequencies and density of states for undoped Hubbard model 
at t = 0.25, U = 2.0, /3 = 20. The results of DMFT and 
the calculation with nonlocal diagram correction (b) are com- 
pared with the reference data obtained for 8x8 lattice QMC 
simulation. 

at Fermi level, without a remarkable anisotropy along 
Fermi surface. At smaller [/, the peak is broadened, with 
maxima near van Hove singularities. The renormalized 
dispersion law -I- ReSi^=o,fc is now also in a qualita- 
tive agreement with numerical data, as the lower panel 
of Figure [5l In these graphs, -|- ReE(j=o,fc is compared 
with the reference data for a 10 x 10 lattice. There is a 
qualitative difference between the results for U = 1.0 and 
U = 2.0: for later case the corrections are quite large so 
that there is a dependence resembling e^^. The superi- 
ority of the result against DMFT should be stressed, as 
there is no fc-dependence of E in the DMFT approach. 

Let us point out the drawbacks of the present results. 
First of all, there is still no perfect quantitative agree- 
ment with the reference data, although the DMFT result 
is improved remarkably. The source of this discrepancy 
becomes clear when the DOS for ?7 = 1.0 is plotted (Fig- 
ure (HI). The pseudogap is much narrower in this case. 
It resembles the situation for U -1-0 at zero tempera- 
ture, then an antiferromagnetic ordering appears due to 
long-range nesting phenomena. Its evident from Figure 
[6] that the calculation with dual diagram (b) does not 
reproduce this pseudogap at all. Back to the results for 
U = 2.0, the pseudogap in our calculation appears to be 
not as deep and not as steep, as it should be (Figure [3]). 
We have tried to take higher diagrams into account and 
found out that it does not help much. We conclude that 
the dual-fermion corrections, as they considered above, 
improve the description of short-range Mott physics, but 
they do not take the long-range antiferromagnetic fluc- 
tuations into account. 

To explain this failure, let us recall the Hubbard model 
with small U at zero temperature. As pointed above, 
our technique passes into weak-coupling diagram expan- 
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FIG. 5: (color online) Momentum dependence for the self- 
energy function at Fermi energy, obtained with diagram (b) 
within the translationally-invariant approximation for the un- 
doped Hubbard model. Data are shown at t = 0.25, /3 = 20, 
for U = 1.0 and U = 2.0. Upper panel: contour plots for 
fc-dependence of the imaginary part of the self energy. Lower 
panel: renormalized dispersion law ek + ReE„=o,fc, compared 
with the reference data obtained for 10 x 10 lattice. The Fig- 
ure has been published previously in the Breief Report [igl. 



sion for (7 — > 0. But it is clear that the weak-couphng 
expansion is suitable for the metallic phase only and 
cannot reproduce the antiferromagnetism, since this is 
a non-perturbativc phcnomenon'^^. Evidently, the dual- 
fermion expansion inherits this property. The best pos- 
sible achievement within this framework would be to ob- 
tain a phase transition, where the corresponding suscep- 
tibility diverges^. 

There are two ways to take antiferromagnetism into 
account. First, one can switch from single-site to clus- 
ter DMFT. Thus the antiferromagnetic phase transi- 
tion maintains the periodical symmetry of the super- 
lattice made of clusters, there is no problem with non- 
analyticity in this case. Indeed, various cluster DMFT 
approachesiiii^ reproduce the antiferromagnetic gap. 
The dual-fermion corrections can be used to improve the 
accuracy of those methods^i^. 

The second option is to stay with the single-site start- 
ing point, but allow for the antiferromagnetic ordering on 
the lattice. In this case, the effective impurity problem is 
spin-polarized. It is known that such an approach indeed 
works quite well already at the DMFT level^^. It can be 
expected, that the dual-fermion technique can effectively 
provide the correction in this case. The next section de- 
scribes such a theory and the corresponding results. 
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FIG. 6: (color online) Density of states for undopped Hub- 
bard model at t = 0.25,(7 = 1.0, /3 = 20. The result of 
the translationally-invariant calculation with diagram (b) is 
compared with the reference data for 8x8 lattice. An anti- 
ferromagnetic pseudogap is pronounced in the reference data 
and does not appear in the approximation. 



B. Undoped case: antiferromagnetic symmetry 
breaking 



For clarity, let us present the explicit expressions for 
this case. The antiferromagnetism means that the primi- 
tive cell is doubled. The dual Green's function, as well as 
other single-electron quantities of the antiferromagnetic 
state, depends on the difference of the two coordinate 
arguments and single spin: G'^^jl j, .,, = G'^j'_y ,. (note 



'dual 



it 



that s' is defined by s and r = j — j'). Given G\ 
is easy to obtain from the formula ([^0)1 . In this 

expression, the spin dependence of S comes from the 
spin polarization of G*'''"'' and, in principle, of the vertex 
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(4) 



However, the numerical result for the latter quan- 



tity appears to be quite noisy. Therefore, we neglected 
the spin polarization of 7^^*^ , performed a averaging over 
spin orientation, and thus operated with the tensor of 
the 'paramagnetic' symmetry. Such a tensor has the two 

independent components 7' = 7isL and 7" = 7s^Ls_s, 
so that the expression (I^O)) becomes 



10 



Y^dual ^ \ / / y~idual j^dual y^dual , ^ \ // // y-idual ^dual y^dual 

uj,r,s c\ 02 / J /cjLJiLj'cj2 iU)Q,uj' uj\ ^r^s^ijj'2^r^s^Lo'^—r^s'^ o2 / j iujujiuj' uj2 iuj2Uj' ujiuj^uji ,r^uj2 ,r,—s^uj',—r,- 



2^2 / J iujuJiuj UJ2 iiJ2UJ iJiuj^ LJi ,r,s^ UJ2 ,r,s^ LU ,—r,s ' p2 

UJ-\-Uj' —UJl -\-UJl UJ-\-Uj' —UJl-\-UJl 




Lattice QMC calculation 

• DMFT 

o diagram (b) 



0.0- 




Lattice QIV1C calculation 




• 


DMFT 
diagram (b) 








-0.2- 


+ 


diagrams {b, 




















-0.4- 






• 


• 

• 




-0.6- 




• 

• 


a 0.00' 




+ 


T J i •■ 

* • 

• 


-0.8- 




• 


^f) -0.05 ■ 




• 


-1.0- 


I 




£ 






-1.2- 


f 




-0.10. 






-1.4- 


• 






• 











0.5 1.0 1.5 



FIG. 7: (color online) Results of DMFT calculation and the 
scheme with diagram (b) , taking the antiferromagnetic order- 
ing into account. The results for local part of the Green's 
function at lowest Matsubara frequency are compared with 
reference data for undoped Hubbard model at t = 0.25, /3 = 
20. QMC calculation at 8 x 8 lattice are used for reference. 



We believe that this approximation is valid, since the 
most important contribution to the symmetry-break 
arises from the spin-polarization of the single electron 
quantities g,A, and E*""^'^ entering the expression for 

The next step is to write explicitly the definition 
^duai = Qdual ~ '^dilai ™ t^i^ momentum space. Here, 
the 2x2 matrices must be used, as the momentum is 
conserved up to Q = (t'",'!'). Let us denote Gf,""5j!^ = 

Gu>,j-j' ,-s)- It is easy to check that the definition 
'^duai = Qdual ~ ^duai '^tays fulfilled with the matrix 

CQdual{a) Qdual(AF) \ 
^dual(AF) ^diial(id) 1 

used for G"*"""', and similarly for Y^duaU Qdual- This gives 
a way to construct G"*""' from a given 1]''"''' and thus 
close the inner iteration loop. The self-consistency con- 
dition remains unchanged, so that the big loop is 
essentially the same. Finally, the exact relationship © 
can be written in the matrix form, giving thus a complete 
description of the antiferromagnetic state. Of course, the 
same treatment with E''"'^' = corresponds to the anti- 
ferromagnetic DMFT. 



FIG. 8; (color online) Imaginary part of the local Green's 
function of undoped Hubbard model at Matsubara frequen- 
cies. The data are shown for J7 = 1, t = 0.25, /? = 20. The 
reference data are compared with the results of approximate 
schemes taking antiferromagnetism into account. The results 
of DMFT calculation, of the scheme with diagram (b), and of 
the approximation taking two diagrams (b), (e) into account 
are shown. Inset shows the deviation of the approximate re- 
sults from reference data. 



Actually, once the antiferromagnetism is taken into ac- 
count, the DMFT result itself already is not too bad. The 
corresponding data are presented in Figure El where we 
show how the Green's function at the lowest Matsubara 
frequency depends on U . At small U, the system is a 
normal Fermi-liquid. There are small corrections due to 
the correlations. Of course, DMFT cannot reproduce the 
anisotropy of the self-energy, but the description of local 
Green's function is pretty good. For large J7, the system 
exhibits a strong antiferromagnetism, which is destroyed 
only at long-range scale. In DMFT, the antiferromag- 
netic ordering appears in this range. The simplest way 
to take the long-range fiuctuations into account within 
DMFT framework is to average over the two antiferro- 
magnet sub-lattices. This eliminates the real part of the 
Green's function. A comparison of IvciG / with lat- 
tice QMC simulations again shows a good agreement (the 
antiferromagnetic regime starts from U « 0.85, as the in- 
set in Fig[7| shows). The largest deviations of the DMFT 
result from the reference data occur in the intermediate 
regime U k, 1. Probably, in this regime the fiuctuations 
are essentially non-local but still mid-range. Therefore 
they cannot be described as a static long-range antifer- 
romagnetic ordering. 
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FIG. 9: (color online) Density of states of the undoped Hub- 
bard model, restored from the data presented in Figure [S] 
The approximate result becomes closer to reference data as 
diagrams (b) and (e) are taken into account. 



The same Fig. [7] presents the result obtained with 
the first nonlocal dual diagram (b). In this calculation, 
we again allow for the antiferromagnetism. The symme- 
try breaking down a almost the same value of U, and 
the magnetization coincides the DMFT result. There is 
however a remarkable correction to ImG^/^ ,,=0- Near 
both limiting cases, the reference dependence is repro- 
duced very well, since the diagram (b) yields a leading- 
order correction to the already good DMFT result. In 
the 'critical' intermediate regime, the situation is not as 
good. However, the correction still behaves regularly and 
shows the correct trend. It is also important that while 
the DMFT data for ImG'7r//3,r=o show a clear kink at the 
transition point, the dual-diagram correction makes the 
curve much smoother. This is certainly more physical, 
because the reference lattice QMC data contains no sin- 
gularities, since there is no true phase-transition. 

We did not found that any particular higher-order dia- 
gram improves the result for G^r/^ ,,=0 significantly. This 
indicates that a large number of higher-order diagrams 
contribute the result. Actually, this is an expectable sit- 
uation near the critical point. However, it was found 
that higher-order ladder corrections give a particularly 
important contribution to the spectral function of the 
system. Let us illustrate this statement, using the data 
for U = 1.0. The Green's function at Matsubara fre- 
quencies for this case are plotted in Figure [51 Since the 
points with dual-diagram corrections are very close to the 
reference ones and can hardly be distinguished, we plot 
also the difference from the reference lattice QMC result 
in the inset of Figure [51 Figure [51 shows the maximum- 
entropy guess for the corresponding DOS. Since the prob- 
lem of analytical continuation of the Green's function to 
the real-frequency axis is known to be ill-posed, we took 
special measures while calculating the density of states. 



The Green's functions are computed with high accuracy, 
and the maximum-entropy analytical continuation is per- 
formed with the same a priory parameters for all curves. 
This ensures that the graphs for the spectral function 
can be compared one with another. The spectral func- 
tion clearly illustrates what is the physical origin of the 
discrepancy between the DMFT and reference data. In- 
deed, since DMFT replaces the nonlocal dynamical anti- 
ferromagnetic correlations with static ordering, it overes- 
timates the antiferromagnetism in the model. Therefore 
the pseudogap appears to be too deep; its shoulders and 
Hubbard bands in the DMFT graph are narrower than 
they should be. The situation is partly improved for the 
diagram (b) : the shoulders and Hubbard bands are closer 
to the reference curve although the estimation at Fermi 
energy looks worse. The serious improvement arises from 
the next diagram of the ladder, as the dash-dot curve in 
Figure [H shows. This is very expectable, because the 
long-range antiferromagnetic fluctuations are exactly de- 
scribed by these ladders. On the other hand it is interest- 
ing to observe from the inset in Figure[9lthat this diagram 
does not improve the result for G.n-/^^r=Oi but makes its 
deviation from the reference data more regular. 



C. Doped Hubbard: Fermi arcs formation and 
flattening of the dispersion law 

Here we present the results obtained with the dual- 
fermion technique for the pseudogap regime, which cor- 
responds to the doping below optimal and relatively high 
temperature. We use the rotationally-invariant approx- 
imation, so the effects of superconductivity and antifer- 
romagnetism were not included in the theory. However 
it turns out that the theory still captures the physics re- 
sponsible for the Fermi arc formation, and yields results 
which compare well to experimental data. 

To make the simulation more realistic we introduce 
the next-neighbor hopping term t' . The parameters of 
the model are U = 4.0, i = 0.25, i' = -0.075, /? = 80. 
The ratio t' /t « —0.3 roughly corresponds to the case 
of YBasCusO,^. The relatively large value of U = 2W 
was taken because there is experimental evidence that 
the system should be a Mott insulator at small doping, 
which requires U > 1.5W ~ 3.0. The temperature used 
roughly corresponds to 100-150 K, which is a proper value 
for the pseudogap phenomena in high-temperature super- 
conducting materials. Most of the results are presented 
on doping level of 14% . 

Figure [TUl presents the results obtained for the self- 
energy S^^fe at the nodal and anti-nodal points of the 
Fermi surface. The position of Fermi surface was defined 
as a maximum of the spectral density. A polynomial ex- 
trapolation for was constructed to obtain the imagi- 
nary part of self-energy at Fermi level. One can observe 
a remarkable difference in the low-energy limit of S„ 
at the nodal and anti-nodal points: the corresponding 
values of ImE^j^o .fc differ approximately by a factor of 
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FIG. 10: Self-energy function of tt' Hubbard model S^^^fc at 
nodal and antinodal points of the Fermi surface at Matsub- 
ara frequencies. Diagram (b) is used for the calculations. 
The data are plotted for 14 % doping tt' Hubbard model 
at t = 0.25, t' = -0.075, U = 4.0, /3 = 80. Upper panel: 
real and imaginary parts of 'Eui.k- Lower panel: ImSuj.k in a 
low-frequency region and its approximation with a 7-th order 
polynomial. 



two. The spectral function Ak — (27r)^^IniGij=o.fc for 
the entire Brillouin zone is mapped in Figure [Til for 14% 
doping. The Fermi surface in the antinodal direction is 
quite diffuse, in accordance with the experimental results. 

It is worth to consider the spatial dispersion of the 
self-energy function. The map of ImS^^o.fc is presented 
in Figure [T^ whereas Fig. [T51 shows the behavior of this 
quantity along the (tt, tt) — (tt, 0) — (0, 0) — (tt, tt) con- 
tour. The data are obtained with a polynomial extrap- 
olation from Matsubara axis. The estimated errorbar of 
the extrapolation procedure is 0.01. An interesting prop- 
erty of the data obtained is that Sc^^o.fc appears to be 
substantially non-local, but still short-range. Actually, 
the data of Figs. [T^ [131 can be approximately described 
by the next-neighbor approximation, that is, the most 
important components of I^uj=a,R. are I]fl=(o.o) , Sfl;=(o.i) 
and £/j=(i_i). The doted line in Fig. [T3| is produced 
with these Fourier-components only, and it is quite con- 
sistent with the initial curve, except the points (tt, 0) and 
(0, tt) where the self-energy is flattened. It is worth to no- 
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FIG. 11: Spectral function Au,=o,k at Fermi level: the calcu- 
lation with diagram (b) and polynomial extrapolation from 
Matsubara frequencies. Parameters of the Hubbard model 
are the same as in Figure [TOl 
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FIG. 12: Imaginary part of the self energy ImE^j^o.fe at Fermi 
level: the calculation with diagram (b) and polynomial ex- 
trapolation from Matsubara frequencies. Parameters of the 
Hubbard model are the same as in Figure 1101 The red line 
indicates Fermi surface. 



tice also that ImEj^^o.fc is maximal there. Interestingly, 
variational cluster calculations^ demonstrate that near 
the nodal point, in contrast with the antinodal one, the 
superconducting gap (that is, anomalous part of the self 
energy) also can be described in the nearest-neighbor ap- 
proximation. 

Figure 1141 shows the effective quasiparticle energy, de- 
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fined by the formula 



4" = Re 



(25) 



The initial dispersion law is shown in the same Figure 
with thin line. One can see an narrowing of the quasi- 
particle band, mainly due to the dTi^^^k/doj term. The 
latter is large due to a closeness to the Mott transition 
point. Another important change is again the flattening 
of the curve near (0, tt) point. 

A flattening of the dispersion curve near the antin- 
odal point was earlier predicted^i^ as due to a non- Fermi- 
liquid behavior when the Fermi energy crosses van Hove 
singularity. The main conclusion is that in the 

strong-interacting regime van Hove point expands to a 
finite region of the Fermi surface, where the dispersion 
law is flattened. The /c-dependence of the self-energy 
and vertex function are of crucial importance for this 
phenomenon. 

Its worth to note that cluster calculation hardly can 
reproduce the result for the van Hove behavior, because 
the flattened region is much smaller that the entire Bril- 
louin zone. 

We also performed calculations for other doping. Fig- 
ure [15] is devoted to ImS at 7% doping. Smaller dop- 
ing makes the system closer to Mott insulator, therefore 
the value of ImE is substantially larger then for the 14% 
doped system (Figs. [TUKTH]). The flattened regions disap- 
pear in this case. However, there is still a clear difference 
between the nodal and antinodal directions in the low 
energy limit: the values of ImE^^^o at these points differ 
by a factor of two. 

Finally, a few words should be said about the region 
near (0, 0) point in Figure 1131 where our polynomial 
fit predicted slightly positive ImE (that corresponds to 
ImG < 0). We argue here that this is merely an artifact 
of the extrapolation procedure. Indeed, as it is discussed 
in section III El negative ImG in our theory could only 
result from a negative residual. However, the graph of S 
at Matsubara frequencies for all /c-points is qualitatively 
similar to whose shown in the upper panel of Figure 1151 
It is obvious these graphs have a negative derivative at 
Fermi energy, so that the residual Z = (i 
be positive. 
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FIG. 13: Imaginary part of the self energy ImE„=o,*: at Fermi 
level: the calculation with diagram (b) and polynomial ex- 
trapolation from Matsubara frequencies. Solid line shows 
the same data as presented in Figure 1121 Dot line is a fit 
with the next-neighbor Fourier components. Arrows mark 
the flattened region at the antinodal direction. Positive sign 
of ImEt^=o,fc in a small region near the (0,0) is probably an 
artifact of the polynomial extrapolation procedure. 
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FIG. 14: Quasiparticle dispersion law, defined from formula 
(|25p (thick line), compared with initial dispersion (thin line). 
Model parameters are the same as in Figures I10I13I Arrows 
mark the flattening of the van Hove singularity. 



IV. CONCLUSIONS 

To summarize, the transformation to dual fermion 
variables completely reconstructs perturbation theory, 
starting with the zeroth-order approximation which is 
accurate in the limits of both very weak and very strong 
interactions. As a result, taking into account just a 
few lower-order diagrams gives quite satisfactory results, 
without having to resumm infinite series of diagrams. 
Starting with DMFT as the best local approximation, 
we are able to take into account nonlocal corrections in a 



regular perturbative way. In contrast with several cluster 
approaches the method is exactly translationally invari- 
ant and allows us to analyze how different parts of the 
reciprocal space are distinctly affected by correlation ef- 
fects. 

This approach can be setup either in phases with long- 
range order (antiferroniagnetism, superconductivity) or 
in phases without long-range order (normal state) by 
not allowing for symmetry breaking. The present ar- 
ticle mostly deals with the latter case. By doing so, 
we could focus on physical effects that are not directly 
related to incipient long-range order. In particular, we 
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FIG. 15: Imaginary part of the self-energy function for the 7% 
doped system. Other parameters of the tt' Hubbard model 
are the same as in Figures 1101 1131 Upper panel: ImE at 
the nodal and antinodal points of the Fermi surface, and its 
polynomial fit at Matsubara frequencies. Lower panel: low- 
energy behavior of ImE at the (tt, tt) — (tt, 0) — (0, 0) — (tt, tt) 
contour. 



showed that the anisotropic destruction of quasiparticles 
and the Fermi surface (at least, as presented in Figs. [TOl 
I14p is not due to precursor effects of antiferromagnetism 
(or superconductivity) as soon as the intermediate and 
strong coupling regimes are entered. Indeed, it is associ- 
ated with quite-short range physics, as illustrated by the 
fact that only the short-range components of the self- 
energy are found to have significant magnitude. This 
observation also provides some support to cluster exten- 
sions of DMFT. 

Although the destruction of quasiparticles in a 
momentum-selective way is adequately captured by this 
approach and associated with short-range correlations, 
more work is required (possibly including symmetry 
breaking and incipient long-range order) in order to reach 
a proper description of the pseudogap formation and of 
its dependence on the doping level and on the t' /t ratio. 

The work was supported by RFFl (grants 08-02-01020, 
08-03-00930, 08-02-91953), DFG (grant 436-RUS-113- 
938-0), FOM (The Netherlands), CNRS and Ecole Poly- 
technique. 



APPENDIX A: EXACT RELATIONS FOR 
HIGH-ORDER CUMULANTS 

Similarly to the exact relationship ^ between the ini- 
tial and dual Green's function, the one-to-one correspon- 
dence between higher-order momenta for the initial and 
dual system can be established. Particularly, the formula 
for the four-order Green's function was presented and dis- 
cussed previously'^^. It was shown that the two-particle 
excitations in the original and dual system are identical. 
Here, we use the generating functional approach, that 
allows us to establish the general structure of such rela- 
tionships for high momenta, and extend the conclusion 
about the two-particle excitations to all collective exci- 
tations, involving an arbitrary number of particles. 

We start from the expression for action ([5]), which in- 
cludes both initial and dual variables. Then we introduce 
the independent variations of initial and dual energy: 

S[c, c*, /, f*;u, v] = S[c, c*J, /*] + Ui2clc2 + V12JU2, 

(Al) 

where u and v are infinitesimal and a summation over 
repeating indices is implied. 

One can see that Taylor series of the functional 



v] —\n e 



^-S[c,c',fJ*:u,v] 



VfVf*VcVc* (A2) 



with powers of u and v correspond, respectively, to the 
cumulants of initial and dual system. We remind that 
the second-order cumulant is the Green's function, and 
higher-order cumulants are proportional to correspond- 
ing vertex parts. For example, the fourth-order cumu- 



lant is 



^1234 - G23G14 + G13G24 {X is the 



two- particle Green's function), whereas the fourth-order 

vertex r[234 = ^^^,622/ Our^i^^' du^, i, ^3'3^4'4- 

To establish a relation between the cumulants, let us 
integrate over /* , / in the previous formula. We obtain 

F[u,v] = Fo[u,v]+ln J e-'^['=''=*;"'"lX>cPc* 



Fo[u,v] = -Indet ||1 + (A - e)a 



S[c, c*; u, v] = S[c, c*] + A^clf^^Cuika + ("12 - M12) C{C2 
M= ((A-e)-i+a-iwa-i)"^ 

(A3) 

Symbol 1 in the second line is the matrix unity, and the 
second term is the product of the corresponding matrices. 
The fourth line reads similarly. 

Last expressions clearly show that the derivatives of 
v\ with respect to u and v are related. A compari- 
son of the first derivatives, for example, allows to repro- 
duce formula The last term of (O comes from the 
differentiation of Fq. 

Let us consider the fourth-order cumulants 



and 



6tl320«41 

First of all we note that neither indices 1 



and 2, not 3 and 4 should coinside, because overvise both 
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cumulants vanish due the Fermi-operator algebra. For 
the case of different indices, the differentiation is quite 
simple and gives, after putting a = 5^^, formula (29) of 
the paper l35l. 



— L1VL2 



dv3'2'dvi> 



-Rz'zRi'i. 



(A4) 



Here L and R are matrix inverse of (A — e)~^g and 
(/(A — e)^"'^, respectively. It should be emphasized that 
this expression does not contain any extra additive terms, 
in contrast to formula ([9]). Formally this is because 



the second derivative 



vanish, as one can check 



straightforwardly. Physically this means that the two- 
particle excitations in the original and dual system are 
identical^. 

There might be also instructive to re-express the last 
formula in terms of vertex function. Putting also a — 
, one obtains 



tI tI j^dual Tj/ r}/ 

1234 — ^ll'^22'^ l'2'3'4'-"3'3-f'4'4i 



(A5) 



where L' = {1 + Y^duaig) ^ and R' = [1 + gY^duai) ^■ 

One can see that the obtained formulas are formally 
valid also for the case of coinsiding indices, when both 
left- and right-hand sides vanish. 

An advantage of the presented approach is that the 
derivation of the formulas for six and higher-order vertex 
parts appears to be literally the same as for the fourth 
order. All the argumentation about the absence of the 
coinciding indices and vanishing of the high derivatives of 
Fq is valid for that case. Therefore formula (|A5|) is valid 
for vertex parts of any order ,just a number of indices and 
multipliers L',R' should be changed. From the physical 
point of view, we conclude that all collective excitations 
of the initial and dual ensemble are the same. 



APPENDIX B: FUNCTIONAL MINIMIZATION, 
RELATION TO DMFT, AND 
SELF-CONSISTENCY CONDITION 

It is clear from the present consideration that a proper 
choice of the hybridization function A is crucial. A 
functional-minimization scheme is suitable to clarify this 
issue. Let us introduce a trial action S[f, /*]. For clarity, 
we put the subscript S at the triangle brackets in this sec- 
tion, to emphasize the the averaging is over the system 
with trial action S. We consider Feynman's variational 
functional 



(Bl) 



A straightforward variation S 
tremum condition 



S + 5S gives an ex- 



<{S~ S)SS >5 = < {S - S) >g< SS>g 



(B2) 



For an arbitrary 5S, this indeed means that the ex- 
tremum of (jBip is delivered by 5 = S*, up to an additive 
constant. In this case (|Bip vanishes. The larger value of 
(|B1[) corresponds to the better approximation. 

There is an important point: since dual action depends 
on A, the condition (|B1[) can be used to determine the 
optimal A. The variation with respect to A gives 



S<S> 



s 



dA 



= 0. 



(B3) 



Here we took into account that variations of S and A 
are independent, so the first two terms of (|Bip do not 
vary with A. As for the last term, it is exactly InZ and 
therefore independent of A as well. 

Now, recalling S[fJ*] = - In / g-^^I^-'^*'/'/'] and sub- 
stituting ([5]), we obtain after certain transformations that 
(|B3[) corresponds to the condition 



g'"'n.u.m 



S' 



(B4) 



Here SgUe is defined by formula (jlOp and Gr=o — 
J2k is local part of the Green's function. While 
deriving these formulas, it is useful to take into account 
that a = g~^ is just a scaling factor standing at /*,/, 
and there is no need to vary this quantity: one can vary 
with respect to A at fixed a and put a = g~^ afterwards. 

Actually, the criterion (jB4| has a very clear meaning: 
local part of the Green's function equals the Green's func- 
tion of the single-site action SsUe, averaged over the fluc- 
tuations of /. Neglecting these fluctuations, one obtains 
just a DMFT condition for hybridization function, that 

is G^,r=0 = gtu- 

To make the consideration more clear, let us flrst 
consider the Gaussian approximation for dual variables, 
S = ^Gdilaif* f ■ Let us show for this Gaussian trial ac- 
tion, the DMFT condition 

GLj,r=0 = guj (B5) 

satisfies (jB4[) exactly (call this statement TI). The proof 
is based on the observation that the condition (jB5p is 
equivalent to the requirement that the local part of dual 
Green's function equals zero. 



-tdual 
'r=0 



0, 



(B6) 



as one can easily check with formulas [T5|) . Further, 
since S is Gaussian, formula (|B6p means that all local 
momenta < /*/, >, < f* f,f* f, >,... equal zero. It 
means that local fluctuations of /, /* are virtually absent, 
therefore < 5™^ [/,/*] >= g^ and jm becomes (iBS]) . 
To obtain a formal proof, one should consider an average 
of the Taylor series for These series starts 

from g^^ whereas the average of any higher term vanishes. 
This profs TI. 
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Next, it is possible also to show that the DMFT 
Green's function is optimal with respect to the varia- 
tions of the Gaussian trial action (call this statement T2). 
With a variation S = -G^lj* f S = -Gj^J* f + 
O/r/2, formula (IB2p becomes 

< is+G^Lirf)fif2 >§=< s+GaLiff >§< nh >s ■ 

(B7) 

The essential point is again that since all local momenta 
of /, /* are vanished because of (|B6p . and the dual po- 
tential V is local in space, all the nonlinearity drops out 
from the (|B7|) . It means that both left- and right- hand 
sides of (|B7|) equal the same value, if —GJuaif*f equals 



the Gaussian part of the dual action. This proofs T2. 

So, we have shown that the DMFT procedure can be 
considered as the Gaussian approximation for the dual 
variables, which is optimal in sense of Fcynman mini- 
mization criterion, with respect to both trial action and 
hybridization function. 

Beyond the Gaussian trial action, an analytical treat- 
ment of the extremal criterion (|Bip is hardly possible. 
Therefore, in the main body of the theory we treat the 
dual system perturbatively, using the diagram series with 
respect to V and the hybridization function defined from 
the condition (fT9|) . 
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